Gate-controlled current switch in graphene 
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Abstract 

We study numerically cross conductances in a four-terminal all-grapliene setup. We show that far 
away from the Dirac point current flows along zigzag directions, giving the possibility to guide the 
current between terminals using a tunable pn-junction. The device operates as a gate-controlled 
current switch, and the electronic properties of graphene are crucial for efficient performance. 
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Graphene, the two-dimensional form of carbon, has been creating a lot of interest not 
only in the physics community but also in the electronic industry, due to its extraordinary 
physical properties These properties include, for example, exceptionally high 

charge mobility, room temperature ballistic transport, ultrahigh thermal conductivity and 
mechanical strength. Provided that the manufacturing questions will be resolved, in the 
long term graphene could be used as material for high-performance nanoelectronic devices, 
even on flexible and transparent substrates. 

Keeping in mind the ultimate goal of building all-graphene circuits, it is important 
to study the electric properties of more complex geometrical structures. Multiterminal 
graphene devices have been studied theoretically in]^, Q]. It was shown, for example, that 
near the Dirac point evanescent modes lead to quantum corrections in the multiterminal 
cross correlations Effects like these are interesting, since they probe the unique elec- 
tronic properties of graphene. One of the most unusual consequences of the quasi-relativistic 
electron dynamics is the Klein tunneling 0, Q, H]. The effect entails that electrons incoming 
to a potential barrier in graphene can couple to the hole states inside the barrier and transmit 
through perfectly, since backscattering is strictly forbidden in clean graphene pn-j unctions 



. Theoretically, graphene nanodevices based on manipulating charge carriers by potential 



barriers have been proposed earlier in 
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In this letter, we study multiterminal conductances in a four-terminal geometry and 
see how the intricate properties of graphene can be exploited even further. We show, for 
example, that the cross conductances between different terminals can be tuned by a step 
potential. In other words, the device performs as a current switch controlled by a top gate. 

The studied geometry is shown in Fig. [TJ The central region (device area) is a rectangular 
graphene island. The size of the the island is given by a pair of numbers {Na, N^), where Na 
(Nz) is the number of outermost atoms on the armchair (zigzag) edge, as shown in the figure. 
The width of the device is W/a = ^/SN^ and the length L/a = 3Na/2^ where a ~ 0.142 nm 
is the carbon-carbon distance. We attach left-right symmetric armchair leads onto the four 
corners of the island and study the cross conductances between the four terminals. One can 
view the first terminal acting as a source of electrons, and our aim is to control the current 
to one of the other leads. 

We calculate the currents numerically in the linear response and zero temperature regime. 
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FIG. 1: (Color online) The setup for an {Na,Nz) device. 



as described by Datta [12|. The retarded Green's function for the device area is given by 



a=l 



(1) 



where Hd is the tight-binding Hamiltonian for the device area and {^a}a=i the self- 
energies describing the four semi-infinite leads. The energy E can be interpreted as the 
overall Fermi level Ep of the system. The tight-binding Hamiltonian of graphene is if = 
— j ^jjc|cj, where tij = t ^ 2.7 eV for nearest neighbors and zero otherwise. The cross 
conductance between leads a and /3 is given by the Landauer transmission formula 
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where s^^ is the scattering matrix describing electron transport from lead [3 to lead a [l3| 



The element of the scattering matrix denotes the probability that a transversal mode 



'a/3 

m in lead [3 is scattered to mode n in lead a, and the trace Tr 
these probabilities. The traces are calculated using the formula 
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FIG. 2: (Color online) The cross conductances as a function of the length Na for a device area with 
width Nz = 52. Fermi level is at Ep = 0.8t. Insets: Schematic representations of the current flow 
for selected aspect ratios (A-C). 

where the coupling matrices are given by the imaginary parts of the self-energies, i.e. 
Ta = i — (Ti^y^. Since we are working in linear response, we refer to "current" and 
"conductance" interchangeably. 

Let us start the analysis by looking at the effects of geometry. The width of device is 
fixed to be A^^ = 52 and we choose leads such that each lead is connected to 13 atoms in the 
zigzag edges. For small Fermi energies, the current is nearly evenly distributed between all 
terminals (see below) . However, for a high Fermi level one starts to see direction-dependent 
effects. 

The results for cross conductances at Ep = 0.8t as a function of the junction length Na 
are shown in Fig. [2l The conductances are given in units of the conductance quantum 
Gq = 2e^ /h. There are nine propagating modes in the leads which is the maximum total 
transmission for leads of this width. Figure [2] shows many interesting features. In a short 
junction, the dominating cross conductance is (ji2, meaning that most of the current travels 
straight from the terminal 1 to the terminal 2 (point A in the figure). Increasing the length 
makes the cross conductance G12 smaller until at Na ~ 50 the cross conductance Gm peaks 
(point B). This can be interpreted as a reflection from the opposite wall, giving the first hint 
that the flow of current is not chaotic by nature. 

When the length of the junction is still increased, the cross conductance G13 starts to 
increase, peaking at Na ~ 100 (point C). Interestingly, this corresponds to an aspect ratio 
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FIG. 3: (Color online) The cross conductances as a function of the Fermi level for a (104,52)- 
junction. The zigzag direction is preferred only near Ep = t. 



of tan~^(M//L) = 30°. This proves that for this high a Fermi level, the current has a strong 
tendency to propagate along the zigzag path from terminal 1 to terminal 3. This is due to 
trigonal warping in the two non-equivalent corner points of the Brillouin zone 

The requirement L = \/?>W means roughly Na = 2A^^, so we fix the geometry to be 
(104,52). The cross conductances as a function of the Fermi level for a (104, 52)-junction 
are shown in Fig. [3l For Ep < 0.5t cross conductances are nearly equal, indicating that 
all terminals are equally probable points for exit. Thus the behaviour of the currents inside 
the rectangular middle area could be described as chaotic. This is not surprising, since the 
relativistic dynamics in rectangular dots is known to be non-integrable [lJ|. Only when 
Ep > 0.5t starts the cross conductance Gis to dominate due to the tendency to propagate 
along the zigzag path. This tendency requires the Fermi level to be rather high, since in 
the Dirac cone all directions are on equal footing, and to get out of this regime the Fermi 
wavelength must be reduced to be of the same order as the lattice constant. 

In the rest of the paper, we study how we can control the direction of the current using 
a potential step. The results for the cross conductances as a function of the height of the 
step are shown in Fig. [H The shape of the potential is V{x) = V\l + exp(— 2x/(i)]^^, with 
a smoothening factor d = a. This is a very steep barrier, and one gets practically the same 
results for (i — i> 0. In Fig. [H one can see three regions of primary interest: (i) For V < 
the current flows mainly to terminal 3 due to the tendency to propagate along the zigzag 
path, (ii) For V ~ 0.8t the right part of the system is at the Dirac point and most of the 
current is reflected from the barrier to terminal 4. The evanescent modes do not contribute 



here, since also the semi-infinite leads in terminals 2 and 3 are at the Dirac point, (iii) For 
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FIG. 4: (Color online) The cross conductances as a function of the step height V for a steep potential 
barrier. The central region is of size (104,52). Insets: (a) Energy diagram in the longitudinal 
direction, (b) Schematic representation of the current flow for V = 2Ef. The barrier works as a 
Veselago lens with refractive index n = — 1, focusing the current to terminal 2. Part of the current 
is reflected to terminal 4. 

V ~ 2Ep the barrier works as a Veselago lens with a refractive index n = — 1, and the 
current transmits to terminal 2. This is the central result of our paper: the current can be 
switched between two terminals using a pn-junction. 

To understand the difference between points (i) and (iii), let us briefly review the principle 
of Veselago lensing ^|. Consider a case in which the left half of the graphene island is at 
chemical potential (jll and the right half at /x/j = — /iL, meaning that there is a symmetric 
np-junction in the middle of the island. Incoming electrons with a propagation angle 6 
can tunnel through the barrier as holes. The propagation angle of the outgoing hole is 
then 9' = —6, since the conservation of transverse momentum requires kcSinO = —kySinO', 
where kc and k^ are the Fermi wave vectors in the n and p regions, respectively. Thus the 
barrier works as a Veselago lens with a refractive index n = — 1, meaning that we can focus 
the charge carriers to terminal 2 instead of terminal 3. Note that the maximum of G12 
is roughly 3/4 of the maximum of G13, indicating an approximate probability of 3/4 for 
transmission through the np-junction. This is in reasonable agreement with the theoretical 
value of T = cos^ (30°) = 3/4. 

Increasing the smoothness of the barrier reduces the cross conductance G12, since 
theoretically, transmission through a symmetric, smooth pn-junction is given by T = 
exp i—TckFdsin^ 9) , where d is the effective length scale of the step and kp is the Fermi 
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wave vector 16!|. But in this case, the current flows to terminal 3 for = and reflects to 
terminal 4 for V" = 2Ef, meaning that the switching effect is not lost but modified. 

It is worth mentioning that we have considered a symmetric structure, in which the 
number of atoms in the leads in the transverse direction is odd. Due to the small size of the 
computational domain, the effects of breaking this symmetry are significant. The precise 
form of the leads is not, however, critical to our results: one can reproduce qualitatively 
similar results in the entire energy range using square lattice leads, when the local potential 
in the leads is chosen to correspond to a good description of contacts |l7l |. 

In conclusion, we have studied cross conductances in a four-probe all-graphene system. 
It was shown that in our setup aspect ratios of taii~^(W/L) = 30° lead to strong turning 
currents due to the preference of current to propagate along a zigzag path, when the Fermi 
level is high. Once this situation is realized, the current can be controlled by introducing 
a potential step in the middle of the island. Smoothening the potential destroys Veselago 
lensing to some degree, but even in this case one can achieve switching between two terminals. 
The fact that the leads were connected to zigzag edges is important, and we leave studies 
on different kinds of setups for future work. 
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